Bending strain-driven modification of surface reconstructions: Au(lll) 
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Strain can affect the morphology of a crystal surface, and cause modifications of its reconstruction 
even when weak, as in the case of mechanical bending. We carried out calculations of strain- 
dependent surface free energy and direct bending simulations demonstrating the change of incom- 
mensurate reconstruction in Au(lll) under strain, in good agreement with recent data. Time- 
dependent strain should cause a sliding of the topmost layer over the second, suggesting an inter- 
esting case of nanofriction. Bending strain could also be used to fine tune the spacing of selectively 
absorbed nanoclusters. 

PACS numbers: 68.35.Bs, 68.35.Gy, 68.35.Md 
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Crystal surfaces often break their ideal bulk-like sym- 
metry, even in the absence of adsorbates, to lower their 
free energy through the so-called reconstructions. Tem- 
perature, adsorbate coverage and (in electrolytes) elec- 
trochemical potential Q cause these reconstructions to 
evolve, either continuously or through two-dimensional 
phase transitions 0- Work on semiconductor surfaces 
j^j-^l indicated that external strain too can drive changes 
in surface reconstruction. In this letter we present theory 
and simulations indicating that reconstruction changes 
at an incommensurate metal surface like Au(lll) can be 
provoked by application of strain, even with the weak 
strains attainable by bending a platelet or slab. When 
compared with very recent data, which appeared while 
this paper was being reviewed, the quantitative agree- 
ment is remarkable. We propose a possible application 
of this fine effect to selective adsorption of nanostructures 
on this surface. Because the strain-induced reconstruc- 
tion changes involve a continuous change of first-layer 
lateral atom density, time-dependent bending strain will 
also drive a kind of incommensurate nanofriction. 

We introduce in a semi-infinite crystal a uniform strain 
\s a f3, A increasing from to 1. The surface free energy 
change, equalling the surface work, is obtained through 
thermodynamic integration: 
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where 7 is the surface free energy per unit area, A the 
area, t™^ 1 the (strain dependent) surface stress, and a, 
(3 span the (x,y) in-plane coordinates. If the surface free 
energies (per unit area) of two different reconstructions 
1 and 2 satisfy 7J < 7°, then in the absence of strain 
state 1 prevails over 2. However, their respective surface 
stress tensors generally differ, and free energies will thus 
vary with different slopes under strain, possibly leading 



to 71 > 72. We are interested in a detectable modifica- 
tion of reconstruction periodicity, particularly if obtain- 
able even with the weak strains provided by mechanical 
bending. 

To explore that, we choose the promising example of 
Au(lll) as our test case. Au(lll) displays a well docu- 
mented (Lxy/3) reconstruction ||, [L » 22.5 in ordinary 
conditions). It is promising because its long incommen- 
surate period L should readily be changeable by a small 
but detectable amount. The (L x -\/3) pattern, commen- 
surate along (112), is incommensurate in the orthogonal 
(110) direction, where the reconstruction consists of a 
succession of of hep and fee stripes, separated by sharp 
soliton-like boundaries, due to an extra row of atoms for 
each period L of reconstruction. X-Ray diffraction exper- 
iments show as temperature is raised from 300 to 800 
K, a continuous shift of the periodicity L from about 22.5 
to 20.9, which confirms the incommensurability. It also 
indicates a denser and denser surface at higher tempera- 
tures, a kind of negative expansion coefficient (probably 
arising, as in the case of Au(100) to compensate for 
a large outwards thermal relaxation of the first layer). 
Here we shall describe how an external strain too can 
change the periodicity L of Au(lll). 

To describe the initial (L x v3) reconstruction we 
use an emp irical many-body classical potential, the glue 
model pQ] , known to describe reasonably well all the 
gold surfaces. In the glue model, the ground state of 
Au(lll) is properly reconstructed, though with a some- 
what shorter T = periodicity L = 11 The differ- 
ence between that and the experimental L = 22.5 is a 
purely numerical discrepancy that can be corrected for, 
and of no fundamental significance. The long range "her- 
ringbone" secondary superstructure also well known on 
this surface |Q is much more delicate involving far smaller 
energies, and will not concern us any further here. All 
calculations and simulations are conducted for thin crys- 
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tal slabs where strain is introduced by simulated bending 
fjJI , as in real-life experiments jj] (although of course our 
theoretical slabs are far thinner). 

Calculations of the surface free energy change under 
strain start with an /V-layer fee (111) slab, its two faces 
(L x \/3) reconstructed with different periodicities L and 
L' . That is obtained by addition of one extra [1 1 2] row 
of surface atoms every L, followed by relaxation 11 1. We 
carry out molecular dynamics simulations of about 1 ns at 
100 K, to obtain by standard methods Jll|] the average 
stress tensor o a f3. From the same simulation we also 
extract a layer-resolved stress profile T a p(zi), defined as 
the force per unit length acting in layer I: 



where i runs over the atoms of layer I, j spans all the 
atoms, and U is the potential energy a function of the 
interatomic distances (this formula is easily modified for 
application to a bent slab |l4|]). The elastic bending work 
is done in our case by r yyi y = [110] being the strained 
direction, orthogonal to the solitons. 
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FIG. 1. Layer-resolved stress profiles of two Au(lll) slabs 
with 16 and 32 layers, subjected to the same bending curva- 
ture. Note the large surface stress and the sub-surface oscil- 
lations. The non linear bulk stress profile is due to anhar- 
monicity. 

The bent slab stress profile obtained for L = 6, L 1 = oo 
(at T — 100 K) is shown in Fig. 1. The smoothly vary- 
ing bulk-like stress in the interior of the slab develops 
at the surface a sharp oscillation, swinging from com- 
pressive in the second layer to tensile in the first layer. 
As shown by comparison of two different slabs with the 
same curvature but different thicknesses, TV — 16, and 32, 
the surface stress of the thinner slab can be immediately 
computed by subtracting off the corresponding bulk-like 



stress of the thicker slab. The bulk stress T yy nlk \z) devi- 
ates from the linearity for increasing strain in a way that 
is very well described by a parabolic fit, the nonlinearity 
due to anharmonicity. The layer-resolved surface stress 
is then obtained as T yy (zi) — T yy (zi) — T yy ulk \zi). Re- 
peating calculations for increasing curvature k we obtain 
the curvature-dependent surface stress. This approach, 
we note, will give free energy differences, fully inclusive 
of the entropic term. Finally the free energy change as a 
function of strain is obtained by integrating |14j : 
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Likewise, the total surface stress can be obtained as 

T a'a s ) — Z/ielayers T af3^^ z ')- 
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FIG. 2. Free energies at 100 K of various (L x y/3) Au(lll) 
surface reconstructions calculated as the strain is changed 
through bending. (The free energy values at zero strain are 
arbitrarily set at their independently obtained T = values) . 
The predicted strain induced drift of L can be extracted from 
the lowest envelope of these curves. 

Fig. 2 displays our main result: the strain-induced free 
energy variation of several different (L x \/3) reconstruc- 
tions of Au(lll) at T = 100 K. The strain dependence 
found for free energy is relatively large, of order 2 meV/A 
(about 1% of the total surface free energy) for 1% strain. 
By comparison with the absolute zero-temperature sur- 
face energy differences at zero strain, independently cal- 
culated (and used in Fig. 2 to fix the unknown relative 
position of curves, thus neglecting an entropy term, hope- 
fully small at low temperatures) we find that 0.18 meV/A 
is roughly the surface free energy difference between pe- 
riodicities L and L + 1 for L around 20 in Au(lll). By 
interpolating Fig. 2 we find that a change of periodic- 
ity from, say, L — 23 to L — 24 can in principle be 
provoked by a compressive surface strain of order 0.2%. 
More generally, strain will cause a periodicity drift de- 
termined by the lower envelope of all the curves: posi- 
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tive (tensile) strain decreasing L, negative (compressive) 
strain increasing it. Eventually for a theoretically very 
large strain L — 00 would be reached, where the recon- 
struction itself disappears. That strain, of order —2% in 
Au(lll) (see below), might however be too large to be 
achieved by bending. 

Our predictions, based on free energies, find a first con- 
firmation in a computer experiment, consisting of a real- 
istic molecular dynamics simulation, based on introduc- 
ing a bending strain in the same slab geometry previ- 
ously used for the surface stress calculation. If we start 
out with two identically reconstructed surfaces, bend- 
ing strain should cause L to spontaneously increase on 
the concave face, and decrease on the convex one. The 
simulation must however overcome a technical problem, 
connected with particle conservation. A change of re- 
construction periodicity L implies a change of first-layer 
lateral surface density, equal to (1 + 1/L) times the den- 
sity of an unreconstructed layer. The initially flat, defect 
free reconstructed surface should, under bending, spon- 
taneously expel atoms from the top layer to form islands 
on the concave face, while opening craters that expose a 
portion of the layer below on the convex face p5| . Such 
phenomena are slow, needing times much longer than 
the typical simulation nanosecond. Because of that, it 
would seem that molecular dynamics simulation should 
be abandoned. A practical way out of this problem |lq ] 
is to include a loose grid of steps in the starting config- 
uration, which amounts to simulate a very close vicinal 
face rather than the original (111) surface. 
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above room temperature, up to about T — 600K . 




FIG. 3. The upper half of the simulation slab, with two or- 
thogonal steps. Periodic boundary conditions require no extra 
steps at the boundaries. The grayscale reflects atomic coor- 
dination: atoms at the solitons between fee and bec regions 
appear darker. 

Our realization, shown in Fig. 3, comprises four sub- 
terraces of heights 0,1,1,2, separated by steps, and con- 
nected one with another across suitably shifted boundary 
conditions Jig). No other steps occur at the boundaries, 
hence the four sub-terraces constitute in reality a single 
large terrace. The terrace size used is 88dx 20\/3d, with d 
the nearest neighbour distance. With this geometry we 
can carry out various simulations, bending an initially 
flat surface and unbending an initially bent surface. In 
order to speed up kinetics, the temperature is raised well 
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FIG. 4. Time evolution of the reconstruction periodicity, 
from a nonequilibrium starting point with L = 16, towards 
equilibrium. Solid line: bent slab with a compressive sur- 
face strain of —4.7%. Dashed line: flat slab with zero surface 
strain. Geometry as in Fig. 3, T= 800 K. Inset: detail of 
the phase difference between the atoms of the first two lay- 
ers, at selected instants of time during the bending. Note 
that the step, indicated schematically at each of the instants, 
has moved slightly to the left in correspondence with the in- 
crease of L (decrease of lateral density). The excess soliton is 
absorbed at the step upper rim between instants (a) and (c). 

Representative simulation results are illustrated in 
Fig. 4, where the continuous drift of reconstruction peri- 
odicity for a surface initially reconstructed with L around 
16 is demonstrated for two different situations, a large 
compressive surface strain, and zero strain respectively, 
the former obtained through simulated bending. The pe- 
riodicity drift from 16 to about 20 and from 16 to about 
11 in the two cases is in good agreement with the free en- 
ergy prediction of Fig. 2 . The inset details the behavior 
of the phase of the topmost atom layer lattice relative to 
the bulk-like second layer (modulo 2ir) as seen on the bent 
slab at selected instants of time. Data are averaged along 
the [1 1 2] direction to reduce the noise due to thermal ex- 
citation: averaging explains why the phase jumps around 
2ir are replaced by a smooth behaviour. We find that the 
periodicity changes are actuated through the appearance 
(or disappearance) of solitons, nucleated preferentially at 
the upper rim of the surface step, where atoms possess 
the lowest coordination. 

Experimentally, attempts to provoke and observe by 
STM surface reconstruction modifications by bending, 
have been carried out recently, precisely on Au(lll). 
Because the strains attained are small they were found 
to affect mostly the secondary structure. An early un- 
published report by Huang et al. described strain- 
induced removal of the wrongly oriented domains of the 
secondary herringbone structure with strains as small as 
1.3 x 10 -5 , with occasional removal of the whole recon- 
struction. Newer data of Schaff et al. jL8| show that the 
secondary structure is again removed but the [L x V3) 
reconstruction survives at least up to a maximum com- 
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pressive strain of 0.4%. By closer examination of their 
data we note that under strain, exerted at an angle of 
30° with respect to the solitons, there is evidence of a 
drift of periodicity L from about 23 to about 25. Wc 
can directly compare that with our calculations. We find 
that the calculated L-dependent surface free energy for 
variable bending is well fitted by a Pokrovskii-Talapov 
form §] AE = -a/L + b/L 3 , where a is the cost of soli- 
tons and b reflects their mutual interaction oc 1/L 2 . Here 
a(e) = a' • (s — s c ) will vanish at a critical strain e c , so 
that L = K(e — Ec)^ 1 ^ 2 , where K = \/3b/a'- Based on 
our computed surface energies and anisotropic surface s 
tress tensor, we derived K and e c for arbitrary strain 
direction, so as to compare the experimental drift with 
calculated periodicity changes (Fig. 5). The agreement 
obtained is very encouraging, fully supporting our quali- 
tative expectations, and demonstrating the applicability 
of the proposed mechanism of surface modification under 
bending strain. 



strain magnitude to actuate the sliding. When the slid- 
ing occurs, a surface nanofrictional contribution should 
appear as a part of the mechanical damping of a Au(lll) 
platelet. Among the applications, there could be the pos- 
sibility to tune through bending strain the spacing be- 
tween deposited nanoclusters. Recently, deposition of 1- 
nitronaphtalene (NN) nanoclusters on a stepped Au(lll) 
surface has revealed JllJ that at the low coverage ad- 
sorption at the step is favored in the fee zones of the 
L x \/?> reconstruction pattern and unfavored in the hep 
zones, and in this way, stripes can be produced. The 
distance among the stripes could in principle be changed 
upon bending, with interesting consequences in the field 
of miniature device fabrication. 

This work was supported through contracts INFM 
PRA NANORUB, and MIUR COFIN. We thank P. Zep- 
penfeld and O. Schaff for bringing their work to our at- 
tention, and for illuminating discussions. 
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FIG. 5. Change of Au(lll) reconstruction period L as a 
function of applied strain. Solid lines, are theoretical predic- 
tions, in the form L = K(e — Ed)" 1 / 2 , for strains orthogo- 
nal (9 = 90°, K = 2.84, e d = -0.0153), oblique (9 = 30°, 
K = 3.34, E d = -0.0210) and parallel (6 = 0°, K = 4.12, 
Ed = —0.0323) to the [1 1 2] oriented solitons . These param- 
eters best fit our simulation results, with a shift in strain, 
to make the L = 23 point correspond to strain e = 0. Data 
points were extracted from the STM images of ref. jl8j , where 
the strain direction was 9 = 30°. 

The strain-induced reconstruction changes described 
above can lead to potentially interesting consequences 
and applications. In a surface that is perturbed by a 
time-dependent strain, there will be as a consequence a 
sliding nanofriction of the incommensurate reconstructed 
top layer of Au(lll) against the regular second layer. 
For example, an oscillatory bending would result accord- 
ing to Fig. 4 in a corresponding variation of L, driving 
an oscillatory sliding of the first layer onto the second. 
A pinned state of the top layer will imply a threshold 
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